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We formulate an exact result, which we refer to as the pumping restriction theorem (PRT). It 
imposes strong restrictions on the currents generated by periodic driving in a generic dissipative 
system with detailed balance. Our theorem unifies previously known results with the new ones 
and provides a universal nonperturbative approach to explore further restrictions on the stochastic 
pump effect in non-adiabatically driven systems. 
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Introduction. The stochastic pump effect (SPE), which 
is the rectification of classical stochastic currents under 
time-periodic excitations, manifests itself in various sys- 
tems For example, it plays a fundamental role in the 
theory of molecular motors [2]. A theory of the SPE, de- 
veloped in the limit of adiabatically slow perturbations 
[s', 'i], provided qualitative insight as well as proposed 
quantitative approaches for pump-current calculations. 

Beyond the adiabatic limit, such a unifying theory is 
yet to emerge. Although many specific systems have been 
studied in great detail, universal results in this domain 
have been very rare. The recently discovered Fluctuation 
Theorems proposed a new approach to the nonadiabatic 
regime [5]. They addressed such questions as work pro- 
duction or probability relations among single trajectories 
j^, but not the magnitude of the SPE. 

In a recent work of Rahav, Horowitz and Jarzynski 
0], a universal result for the non-adiabatic regime, called 
the "no-pumping" theorem, was formulated. One can 
parametrize kinetic rates of an arbitrary Markov chain 
with detailed balance conditions, so that for any pair of 
sites i and j these rates are written as kij = ke-^^~^^\ 
where Ei can be called the depth of a potential well z, and 
Wij = Wji is called the size of the potential barrier i-j. 
Energy scale is ksT = 1. The theorem [7| states that to 
generate rectified currents during a cyclic process, both 
well depths and barrier sizes must be varied. 

In this Letter, we derive generic restrictions on the 
SPE, that include previously found theorems [1, as 
a special case. We show that there is a wide class of 
stochastic models on graphs, where nonadiabatic but pe- 
riodic modulation of some parameters leads to a zero 
time-averaged flux through any link. We also predict re- 
strictions on the values of nonzero pump currents when 
they are allowed. Such restrictions are sensitive to the 
topology of the graph representing the system. Thus, we 
claim that it is possible to determine whether two given 
links on a graph belong to the same loop by merely mea- 




FIG. 1: Five-state graph with one loop [2342]. 



suring the presence of the SPE. 

Representative examples. Consider a particle moving 
according to Markov chain rules on a graph with a fi- 
nite number of states and the detailed balance. First, 
we note a trivial observation, that periodic variation of 
parameters on a tree-like graph would not produce a net 
current because any flux that passes through any given 
link would eventually return through the same link after 
driven parameters return to the initial values. Pump cur- 
rents are possible only on graphs with loops, such as the 
one shown in Fig. [H Note, however, that the previous 
analysis still applies to some of the links. Namely, the 
integrated over time fluxes through links 1-2 and 4-5 in 
Fig. [T] must be zero. 

A less trivial problem is whether there are general 
conditions (beyond the situation discussed in [7]) under 
which the flux through any link belonging to the loop in 
Fig [H such as the link 2-3 can be zero. For the model 
in Fig [T] we claim, that if one varies the rates related 
to the links 1-2 and 4-5 cyclically but otherwise arbi- 
trarily, while other rates remain constant and satisfy the 
detailed balance condition, then the time integrated cur- 
rent through any link on the loop will be zero. To prove 
it, we split the set of links in Fig. [T] into two subsets {1- 
2, 4-5} e Xi and {2-3, 3-4, 2-4} G Xq. Consider the 
evolution of probabilities pi on three sites connected by 
the links in Xq, i.e. the sites with indexes i = 2,3,4. 
Conservation laws require that 



P2 = jl2(t) -j23{t) -j24{t), 



(1) 
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FIG. 2: The total flux passed through the link 2-3 by time t. 
Variable parameters are Ei(t) = ^i+^sin(a;t), and W45(t) = 
V^45 + Vl/cos(u;t), for t e (0,27r/u;), and Ei{t) = Ei, W45{t) = 
W45 + for t > 27r/cj, and we distinguish between Ei{t) and 
El etc. Choice of constant parameters isa; = l, ^=1, VK = 
2, El = -1, E2 = 0.2, E3 = -0.05, E4 = 0.1, Eb = -0.25, 
W45 = 0.3, W12 = 0, W23 = 0.2, 14^34 = -0.2, 1^24 = 0.25, 
Wij = Wji for any i-j. 

where jij{t) is the current passing through the link i-j. 
Requiring the final probability distribution to be equal 
to the initial one, we reach the conditions 

/ dtpi= f dtji2{t)= f dtj45(t)=0, (2) 

Jo Jo Jo 

where the upper limit of the integration depends on the 
choice of the driving protocol. In case of a single localized 
pulse that drives kinetic rates on links in Xi, T must be 
formally infinite in order to allow the system to relax to 
the equilibrium. In case of a steady periodic driving with 
the period r, if the system already reached the steady 
regime with pi{t) = pi{t -\- r)^ one can choose the upper 
integration limit T = r. 

The current through any link can be formally written 
in terms of instantaneous site probabilities and the ki- 
netic rates, jij = kjiPi — kijPj. Denote pi = pi{t)dt. 
Integrating ([T|) over time, using conditions ([2]) one arrives 
at a set of equations, 

-(^2 + k42)p2 + k23p3 + ^24p4 = 0, 

-{k23 + ^43)^3 + ^2^2 + ks4P4 = 0, (3) 

-(/C24 + k34)p4 + k42p2 + ^^^43^3 = 0. 

The set of equations ([3]) for pi coincides with the one for 
probabilities on a 3-state Markov chain with the same 
topology and rates as for a subset Xq at equilibrium state. 
Since all rates on this subset are time-independent and 
satisfy the detailed balance, the solution of ^ is pi = 
Ce~^% i = 2,3,4, where C is a constant that depends 
on the details of the driving protocol, but is equal for all 
three states on the loop. The total flux passed through 
e.g. the link 2-3 then reads 

J23(T) = / j23{t)dt = C(fc32e-^^ - fc23e-^^) = 0. (4) 
Jo 



^23 

OS 
0.06 








— ^Ei(t),w34m 

-*-Ei(t),W24(t) 




0.04 


■ j\ ^ 
if \ ^ 






0.02 


■J \ 







-0.02 


'& O 



5 10 15 20 25 30 t 



FIG. 3: The total flux passed through the link 2-3 by time t. 
Wij is varied along one of the loop links. Variable parameters, 
indicated on a graph legend, change with time according to 
E^{t) = E^ + Esm{oot), and Wij{t) = W^J +W cos{oot), for t G 
(0,27r/cj), and Ei{t) = Ei, Wij{t) = Wij + Vl^ for t > 27t/lv. 
Constant parameters are as in Fig. [2] 



This concludes our proof that varying the kinetic rates 
outside the loop only, does not lead to a net time averaged 
flux through any link on the graph in Fig. [TJ 

Numerical check. To test our predictions we performed 
numerical simulations for the graph shown in Fig. [H 
In all examples we assume the kinetic rates to satisfy 
the detailed balance condition with the parametrization: 
kij = kex-pdEj — Wij]/kBT)^ [7]. We chose the energy 
scale so that /c = 1, ksT = 1. Fig. [2] shows the flux 
passed through the link 2-3 on the loop, when only pa- 
rameters El and W45 outside the loop are changing with 
time. After completion of the external driving, the sys- 
tem relaxes with a zero net current through the link 2-3, 
as we predicted. 

Our numerical results show that if barriers are var- 
ied on the loop but potential depths Ej are varied only 
on external sites then nonzero pump current through 
the loop links appears, as shown in Fig. [3l However, 
the opposite is not true. Namely, if one varies well 
depths on the loop sites but barriers are varied along 
external links, numerically we always found no overall 
pump flux, as we show in Fig [H The latter result is 
also easy to understand by introducing the quantities 
fi = dt[e^^^^^pi{t)], i = 2,3,4. The same analysis 
as before leads to the equilibrium master equation for fi 
with equal forward/backward rates connecting any pair 
of sites, and a constant solution, ensuring the zero time- 
averaged current through any link on the loop. The pre- 
viously found theorem [7] can be proved for any graph 
by the same steps and change of variables. We are now 
in a position to frame our result in a form of a general 
and rigorous mathematical statement. 

Pumping- Restriction Theorem (PRT). For a graph X 
we denote the vector spaces of time averaged popula- 
tions p G Cq{X) and time-averaged currents J G Ci(X), 
so that p = {pa}, J = {Jab} with Jba = -Jab and 
Jab — when the nodes a and h are not connected by 
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FIG. 4: The total flux passed through the link 2-3 by time t, 
when only Ei were allowed to vary on one of the loop sites 
and Wij is varied along external links. Variable parameters, 
indicated on a graph legend, change with time according to 
Ei{t) = Ei-\-Esm{ujt), and Wij(t) = W^J -\-W cos{ujt), for t G 
(0,27r/u;), and Ei{t) = Ei, Wij{t) = Wij + Vl^ for t > 27r/cj. 
Constant parameters are as in Fig. [2l 

an edge. We further introduce the boundary operator 
d : Ci{X) ^ Co{X) by {dJ)a = EbJba- Let Hi{X) be 
the subspace of physical time- aver aged currents, i.e. sat- 
isfying the continuity condition, and Hq{X) represent the 
space of populations that would be constant within the 
connected components of X. Note that the conjugate op- 
erator 9''' : Cq{X) Ci{X) has a form {d^p)ab = Pa — pb- 
In detailed balance the rates are given by kab = Qab^^^ ^ 
where Qab = 9ba = A:e~^^^ is referred to as a metric on 
X. The Euler theorem claims ^] 

dimCi — dimCo = dimi^i — dimi^o- (5) 

Consider a partition X = Xq U Xi, where Xq consists 
of edges (and adjacent nodes) where the rates are given 
by kab{t) = gab^^^^^\ whereas Xi represents the rest of 
the edges with arbitrary rates kabif)- Note that Xq flXi 
does not have any edges, whereas its nodes provide the 
currents between Xq and Xi. The PRT claims that 

(i) for the described periodic driving the generated 
pumped current is restricted to a vector subspace V such 
that J G F C i^i(X), with the dimension 

dimV = dimi^i(X) - dimi^i(Xo) = dimCi(Xi) - 

-dimi^o(Xo) - dimCo(X) + dimCo(Xo) + 1. (6) 

(ii) given a set of links with driven barriers, the embed- 
ding V C Hi{X) is totally determined by the (time- 
independent) metric Qab on Xq. 

The statement (ii) implies that if there are constrain- 
ing equations that determine relations among possible 
pumped currents through different links on Xq, then co- 
efficients in these equations will depend only on the met- 
ric on Xq. Note also that the choice Xq = X repro- 
duces the second no-pumping theorem of [7]. The PRT, 
however, claims more: starting with the no-pumping sit- 
uation with at least one of driven, and also driving 




FIG. 5: A six state Markov chain. 



the barriers at a certain number n of links, the num- 
ber of independent generated currents may not exceed n, 
i.e. dimV < dimCi(Xi), that follows from second equal- 
ity in and obvious inequalities dimi^o(Xo) > 1, and 
[-dimCo(X) + dimCo(Xo)] < 0. Each of the driven links 
can be viewed as either responsible for an independent cy- 
cle or for connecting two disconnected parts. Therefore, 
the PRT can be interpreted as the claim of the number of 
independent generated currents to be equal to the maxi- 
mum number of driven barriers, which removal does not 
split the remaining graph into disjoined components. 

The proof of PRT is based on the Master Equation and 
the expression for the current 

m = dm, m=0e^^')p{t), (?) 

where the second equality is valid on Xq. Representing 
the current j{t) — jo{t) -\- ji{t) as the sum of the Xq and 
Xi components and averaging (integrating) Eq. ([7j) over 
time we obtain for the time-integrated quantities 

dgd^f = -C. C = dJi\xonx,. J^=0f. (8) 

where / = e^^^^p{t), Ji/o = Jq ji/o{t) and C ^ 
i^o(Xo n Xi) describes the flux passed between Xq and 
Xi. Since the operator dgd^ can be viewed as the dis- 
crete Laplacian in Xq associated with the metric ^, the 
flrst two equations in ([8|) have a solution if and only if the 
total flux entering any connected component of Xq and 
Xi is zero. Let Vo C Ho{Xo HXi) is the vector subspace 
of that satisfles these conditions. The solution of the 
flrst equation in ([8]) is unique up to an additive constant 
distribution, which does not affect the value of Jo, i.e., 
the latter is uniquely determined by ^ G Vb and metric 
on Xq. On the other hand, given ^ G Vb, the component 
Ji is deflned up to a current that circulates completely 
within Xi and is represented by an element of Hi{Xi). 
Thus we proved that the metric on Xq and topology of 
Xi determine the restrictions on J, which implies the 
statement (ii) of PRT. In addition our analysis implies 

dimV = dimFo + dimi^i (Xi ) . (9) 

To derive the explicit expression (|6|), which completes 
the proof we combine Eqs. ([9]) with the identity 

dimi^i(X) = dimi^i(Xo) + dimi^i(Xi) + dimVb, (10) 

which results in the flrst equality in (j6|). The second 
equality follows from ([5|) and that dimi^o(X) = 1. Note 
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that the identity (p!Q|) has a very simple physical meaning: 
For any physical current J G Hi{X) on our graph we can 
identify the current J^^^^ G Vq that flows from Xq to Xi, 
and once the exchange current J^^^^ is identified, the 
complete current J is defined up to the currents J^^^ G 
Hi{Xq) and J^^^ G Hi{Xi) that circulate strictly within 
and Xi, respectively. 

We further illustrate the PRT using a graph in Fig. O 
with dimCi(X) = 7 (the number of links), dimCo(X) = 6 
(the number of sites), and two independent loops, e.g., 
[2342] and [2462] that form a basis in Hi{X) with 
dimi^i(X) = 2. If only one barrier is driven we can 
generate not more than a ID subspace of currents, i.e. 
dim]/ < 1. If the barrier at links 1-2 or 4-5 are driven, we 
are at the no-pump situation, since d\mV = (the sub- 
graph Xo obtained upon elimination of the link 1-2, or 
4-5, has two loops, hence dimi^i(Xo) = 2). Driving the 
barrier at any other single link yields dimV = 1, since the 
graph Xi in this case is connected. When a pair of barri- 
ers is driven we have dimV = for {1-2, 4-5}; dim]/ = 1 
for {2-3, 3-4}, {2-6, 6-4}, and when one of 1-2, 4-5 and 
one of the rest are driven. In all other cases we have 
dim]/" = 2 (all currents may be generated). Consider 
in more detail the {2-3, 3-4} driving. The subgraph Xi 
includes these two links and the vertices {2, 3, 4}, Xq in- 
cludes the rest of the links and the vertices {1, 2, 6, 4, 5}, 
whereas Xq Pi Xi has the vertices {2,4} and no links. 
The exchange goes through two nodes {2,4}, the ex- 
change currents satisfy relations J23 + J43 = that yields 
dim]/b = 2 — 1 = 1, i.e. while allowed elements of XoflXi 
have the form {p2, P4}, elements in Vq are restricted to be 
of the form {J, — J}. The subgraph Xi has no loops, and, 
therefore, no internal currents, which yields dimV = 1, 
which agrees with the PRT. 

We also note that similar arguments lead to continu- 
ous (Langevin dynamics) counterparts of the no-pumping 
theorems of [7]. Let g^^ (x) be a metric in a com- 
pact oriented manifold M that describes bath-induced 
fluctuations/dissipation. If parametrized as g'^^{x) = 
h^^{x)e^^^''^\ where h is time-independent and ]/ is a 
periodic in time potential, then {x) = 0. 

Conclusion. Many emerging mesoscopic devices, in- 
cluding molecular motors [3] and nanoscale electronic 
circuits [9] can be modeled as discrete connected enti- 
ties with stochastic transitions among different states. 
The control over such new devices is impossible with- 
out deep understanding of non-adiabatic strongly driven 
regimes of their operation. Our result is a step toward 
such a theory. The PRT, presented here, determines the 
restrictions in the space of pumped current values on a 
graph. No-pumping conditions follow as its special con- 
sequences. The restrictions on the pump current space 
suggest, for example, that an application of a periodic 
stimulus can be used to induce a localized rectified effect 
without perturbing the whole circuit on average even if 



all its components are connected. It should be useful 
to explore how the PRT is modified by quantum effects. 
Violation of the "no-pumping" conditions due to quan- 
tum corrections can be employed to detect the quantum 
pump effect [10] by cooling the electronic circuit in the 
"no-pumping" regime down to the quantum domain. 

Another possible application for this work is the re- 
construction of the stochastic network topologies, e.g. 
in biochemical reactions. Standard measurement tech- 
niques, such as those based on the linear response, appear 
insufficient. One would expect to find a nonzero signal 
anywhere on an ergodic Markov chain in response to a 
time-dependent current inducing perturbation. However, 
we showed that measuring rectified current in response 
to external periodic stimulus can help one to identify 
whether or not two given links belong to the same loop. 
One of the perturbed links must belong to a common 
loop with the measured one in order to observe the SPE. 
Even without a quantitative understanding of the data, 
detecting only the presence of the SPE can be sufficient 
to deduce the topological structure of the network com- 
pleting only a small number of measurements. 
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